Investigating polaron transitions with polar molecules 
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Abstract 

We determine the phase diagram of a polaron model with mixed breathing-mode and Su- 
SchriefFer-Heeger couplings and show that it has two sharp transitions, in contrast to pure models 
which exhibit one (for Su-Schrieffer-Hccgcr coupling) or no (for breathing-mode coupling) transi- 
tion. Our results indicate that the physics of realistic mixed polaron models is much richer than 
that of simplified models. We then show that ultracold molecules trapped in optical lattices can 
be used to study precisely this mixed Hamiltonian, and that the relative contributions of the two 
couplings can be tuned with external electric fields. The parameters of current experimental set- 
ups place them in the region where one of the transitions occurs. We propose a scheme to measure 
the polaron dispersion using stimulated Raman spectroscopy. 



Introduction: Polarons, which are the low-energy dressed quasiparticles appearing in the 
spectrum of particles coupled to bosonic fields, have been of broad interest in physics ever 
since their first study by Landau [1]. The generic Hamiltonian for a single polaron problem: 



includes the kinetic energy of the bare particle (first term) , the energy of the bosonic mode 
(second term) and their interaction (third term), where the particle scatters by absorbing 
or emitting bosons. 

There are two mechanisms for the particle-boson coupling, as the presence of bosons can 
change (i) the potential or (ii) the kinetic energy of the bare particle. The former always 
leads to a vertex gg independent of the particle's momentum k, while in the latter case 
the vertex gk,q depends explicitly on both momenta. As an example, consider interactions 
between electrons and phonons. Vibrations of nearby atoms modulate the potential energy 
of an electron. Much studied examples of such type (i) interactions are the Holstein [2] and 
Frohlich [3] models. At the same time, by modulating the distance between lattice sites, 
lattice vibrations also affect the hopping integrals. Such effects are described by type (ii) 
models like the Su-Schrieffer-Heeger (SSH) model [4]. 

Most of the early polaron studies focused on type (i) models, in particular on the search 
for a self-trapping transition whereby at strong coupling the bosons create a potential well 
so deep that it traps the polaron. In the absence of impurities, however, this idea has proved 
to be wrong: type (i) polarons have a finite effective mass for any finite coupling strength. 
The absence of transitions for type (i) models has been demonstrated analytically [5] . These 
models have a smooth crossover from light, highly mobile polarons at weak coupling to 
heavy, small polarons at strong coupling. 

This standard view of the polaron as a quasiparticle that becomes heavier with increased 
coupling is now strongly challenged by results for type (ii) models. Recent work has shown 
that in such models, the polaron can be lighter than the bare particle, since the bosons 
affect the particle's hopping so it may move more easily in their presence [6-8] . The boson- 
mediated dispersion can be quite different from that of the bare particle, for example it 
may have a ground state with a different momentum [7]. If this happens, a sharp transition 
should occur when the boson-mediated contribution to the dispersion becomes dominant. 
Indeed, such a transition was recently predicted for the single SSH polaron [8, 10]. 
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These results are interesting not just for challenging a long-estabhshed paradigm, but 
also for raising the question of what happens in realistic systems, where a mix of both 
types of coupling is generally expected. For example, what is the fate of the transition 
in a mixed model if the coupling is varied smoothly from gk^q to Qq'! This also makes it 
highly desirable to find systems described by such mixed Hamiltonians but where, unlike in 
solid-state systems, the various parameters can be tuned continuously so that the resulting 
polaron behaviour can be systematically investigated. 

In this Letter we elucidate the evolution of this transition as the coupling interpolates 
between type (ii) SSH and type (i) breathing mode (BM) [11]. Surprisingly, we find that 
the phase diagram has two sharp transitions, and that these may occur even when the type 
(i) coupling is dominant. This shows that polaron physics is much richer than generally 
assumed. We then show that this mixed Hamiltonian describes polar molecules trapped in 
an optical lattice, and moreover, that the parameters of current experimental set-ups place 
them in the region where a transition is expected to occur. Thus, experimental confirmation 
of these transitions is within reach. Furthermore, we propose a detection scheme equivalent 
to Angle- Resolved Photoemission Spectroscopy (ARPES) in solid-state systems [12], which 
directly measures the polaron dispersion and can therefore pinpoint the transition. Finally, 
as we discuss in conclusion, such experiments hold the promise not only to clarify many 
other effects of interplay between type (i) and type (ii) couplings on single polarons, but 
also to investigate the finite concentration regime to look for quantum phase transitions. 

Phase diagram: Consider the particle-boson coupling 

2i 

9k,q = -y= {a [sin(A; + - sin(A;)] + ^ sin(?)} , (2) 
V A 

where A^ is the number of lattice sites in a one-dimensional chain. The (/c, g)-dependent part, 
with energy scale a, describes SSH coupling [8] while the /c-independent part, with energy 
scale ^, describes BM type modulations of the site energy by bosons [11]. We assume 
Einstein bosons = Q and a free particle dispersion = +2tcos{k) with t > [13]. 
Following Ref. [8], we define the effective SSH coupling strength A = 2a^/{thQ) and the 
adiabaticity ratio A = hfl/t. In addition, we introduce R = (3/a to characterize the relative 
strength of the two types of coupling. 

The SSH polaron {R — 0) was predicted to undergo a sharp transition at a value A* 
[8]. Its physical origin is simple to understand in the anti-adiabatic limit ^ ^ 1 where 
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FIG. 1: (color online). Phase diagram A vs. A at fixed R. Symbols show the location of the 
transition if i? = (blue circles) and R = —0.5 (red circles). The three shaded regions represent 
molecules with large (a, LiCs), intermediate (b, RbCs) and small (c, KRb) dipoles, and dressing 
schemes 1 {R = —0.5 for a) and 2 (i? = for b and c), respectively. For each region, the three 
lines correspond to lattice constants of 256 nm (lower bound), 532 nm (dashed line) and 775 nm 
(upper bound, not shown for region (c) because it is off-scale). For each curve varies between 
1 and 100 kHz. 

the SSH coupling leads to an effective next-nearest neighbor hopping i ^ z + 2 of the 
particle, by first creating and then removing a boson at site i + 1 [8]. Its amplitude is 
t2 = —a'^/{hQ) = —At/2 < 0, so its contribution — 2t2Cos(2A;) to the total dispersion has 
a minimum at 7r/2, unlike the bare dispersion which has a minimum at vr. If 4|t2| = t, 
corresponding to A* = ^ in the limit of A ^ 1, a sharp transition marks the switch from a 
non-degenerate ground state with momentum kgs = tt (for A < A*) to a doubly-degenerate 
one with \kgs\ — )■ | (for A > A*). As A decreases the number of phonons in the polaron 
cloud increases. This renormalizes both hoppings t — )■ t*, ^2 — ^ ^2; so A* changes smoothly 
with A as shown by the blue circles in Fig. 1 (see also Fig. 4 of Ref. [8]). We also plot A* 
for R = —0.5 (red circles), showing that the sharp polaron transition persists for coexisting 
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FIG. 2: (color online). Phase diagram A vs. R at fixed A, showing two sharp transitions: one 
from a non-degenerate ground state with kgs = vr to a doubly-degenerate ground state with < 
l^gsl < T^, and the second back to a non-degenerate ground state with kgs = 0. The results are for 
A = hQ/t — )• oo (red lines) and for A = 5 (blue circles). 

type (i) and type (ii) couplings. These results were generated with the Momentum Average 
(MA) approximation, specifically its variational flavor where the polaron cloud is allowed 
to extend over any three consecutive sites [7, 8]. For A > 0.3, MA was shown to be very 
accurate for both SSH and BM couplings [8, 11]. 

To understand the evolution of A* with R, consider again the limit 1. In addition to 
the second-nearest neighbor hopping t2, there is now also a dynamically generated nearest- 
neighbor hopping ti = 2aP/{hQ) = RXt. This describes processes where the particle hops 
from site i to i + 1 leaving behind a boson at i (SSH coupling) followed by absorption of 
the boson while the particle stays at i + 1 (BM coupling); or vice versa, hence the factor 
of 2. The total nearest-neighbor hopping is thus t* = t — ti, and the transition now occurs 
when A\t2\ = t* -> A* = l/{2 + R). Thus, for i? < 0, the interference between the SSH 
and the BM couplings results in a larger effective t* leading to a larger A*. In particular, 
A* — )■ oo for R < —2, i.e. no transition occurs here. The lack of a transition is not surprising 
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when it! — > — oo, since here the BM couphng is dominant and pure Qq models do not have 
transitions [5] . Our results show that for mixed SSH+BM coupling, the switch from having 
to not having a transition occurs abruptly R — —2 \i A ^ 1. This value must change 
continuously with A, therefore we expect this switch to always occur at a finite R. 

This is confirmed in Fig. 2, where we plot A* vs. R for A ^ oo and ^4 = 5. Surprisingly, 
we find not just the transition at A* ~ + 2), but also a second one which marks the 
crossing to a ground state with kgg — 0. Its origin is also easy to understand in the anti- 
adiabatic limit: if R\ > 1, T is negative and favors a ground state at k^^ — instead of 
/Cgs — TT. For A^ 1 this second transition is at A* = l/(i? — 2) if i? > 2. At finite A, it 
moves towards smaller (i?, A) values, sec Fig. 2. 

Interestingly, this shows that for R — )■ +oo there are two nearby transitions for the shift 
kgs = TT to /cgs = 0. This seems to contradict the proof that a type-(i) Hamiltonian cannot 
have transitions [5], however, even for — > oo this is a mixed Hamiltonian if A 7^ 0. The 
transition is indeed absent if a = 0. This is an example of the rare occurrence where a 
perturbatively small term has a large effect on the behavior of the system. 

Cold molecule implementation: Polar molecules in optical lattices can be used to imple- 
ment Hamiltonian (1) in a wide region of the parameter space. Specifically, we consider 
molecules prepared in the ro- vibrational ground state of the spinless electronic state 
and trapped on an optical lattice in the Mott insulator phase, as recently demonstrated 
experimentally [14, 15]. We assume that there is at most one molecule per lattice site. 

The dipole-dipole interaction between molecules in different sites can be modified by 
applying a DC electric field E = Edqt, [16-21]. Here we consider two schemes for dressing 
the rotational states of molecules with electric fields that are relevant for polaron observation, 
scheme 1 involving a DC electric field only, and scheme 2 involving combined optical and 
DC electric fields. For the former, we define the two-state subspace \g) — |0, 0) and |e) = 
|i,0), where |7V,Mjv) denotes the field-dressed state that correlates adiabatically with the 
field-free rotational state \N,Mn)- N is the rotational angular momentum and Mjv is 
the projection of N along the electric field vector. In this basis we define the pseudospin 
operator c| = \ei){gi\ that creates a rotational excitation at site i. This excitation (the "bare 
particle") can be transferred between molecules in different lattice sites with an amplitude 
tij — jUij (1 — 3cos^ O), where Uij — d? /\vi — rjp, d is the permanent dipole moment, is 
the position of molecule in site i, © is the polar angle of the intermolecular separation vector. 
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7 = IJ'lg/(P < 1 is the dimensionless transition dipole moment that depends on the strength 
of the DC electric field. The excitation hopping amphtude is finite even for vanishing field 
strengths. The field-induced dipole-dipole interaction shifts the energy of the state |ej) by 
Di = ^jDij. Here Dij = —nUij (1 — 3cos^0), where k = \^g{^e — fJ'g)\/d'^ and /ig(/ie) is 
the induced dipole of the ground(excited) state. Dipolar couplings outside this two-level 
subspace are suppressed when the electric field separates state |e) from other excited states. 

The free quasiparticle dispersion is ek — Sq + 2tcos{k) where the site energy is Eq — 
hbUeg + DQ, with the single-molecule rotational excitation energy huUeg ~ 10 GHz and t = ti2- 
The center-of-mass vibration of molecules in the optical lattice potential is coupled to their 
internal rotation through the dependence of Uij on — r^. For harmonic vibrations with 
linear coupling between internal and external degrees of freedom [18, 22], the boson term in 
Eq. (1) describes lattice phonons whose spectrum depends on the trapping laser intensity 
and the DC electric field [22]. Here we consider weak DC fields and moderate trapping 
frequencies which give a gapped and nearly dispersionless (Einstein) phonon spectrum with 
frequency Q. 

With these definitions, the system is described by SSH and BM-like couplings with the 
energy scales a = —3{ti2/ci.L)\/~^J^^rnIl and /3 = —3{Di2/aL)\/^/'2mQ, respectively, where 
m is the mass of the molecule and is the lattice constant. The ratio R — 13/a — 
Hg{He ~ l^g)/l^% is independent of the intensity of the trapping laser or of the orientation of 
the array with respect to the DC field. In the field dressing scheme 1 we have \R\ < 1/2 
for dEDc/Be < 1. In the combined AC-DC dressing scheme 2, the same DC field strength 
and orientation is used as above, but an additional two-color Raman coupling redefines the 
two-level subspace as \g) = ^/a |0,0) + -\/l — a |2,0) and |e) = Vb [1,0) + y/1 — b |3,0) 
(details in the Supplementary Information). This dressing scheme can be used to effectively 
enhance the hopping amplitude by a factor / > 1 yielding a ^ fa, without changing the 
value of nor the phonon dispersion. When using this dressing scheme, any point in the 
phase diagram transforms as A — > /A and A — > A/f, thus shifting the system towards 
stronger SSH couplings. 

The frequency of lattice phonons in a ID array is f2 = (2//i) VVo-Er where Vq is the lattice 
depth and = h^n'^ /2ma\ is the recoil energy. The particle-boson coupling can thus be 
written as A = 18(£'ij/i)(7rA)~^. The shaded regions in Fig. 1 show accessible points in the 
polaron phase diagram (A, ^4) for LiCs, RbCs and KRb molecules, illustrating the fiexibihty 
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in varying the Hamiltonian parameters when using the two field dressing scenarios and 
different experimental settings. Figure 1 shows that the transition characterized by the shift 
from a non-degenerate ground state kgg = tt to a degenerate ground state < \kgs\ < tt 
can be studied using molecular species with moderate dipole moments such as RbCs, in 
lattices with a site separation ~ 500 nm. However, the transition is easier to observe for 
molecules with large dipole moments such as LiRb and LiCs. For weakly dipolar molecules 
such as KRb, the strong coupling region can be achieved using ql < 500 nm. 

The most direct way to detect the transition is to measure the polaron dispersion. We 
propose the stimulated Raman spectroscopic scheme illustrated in Fig. 3 to achieve this goal. 
We consider a one-dimensional array initially prepared in the absolute ground state l^f) = 
. . . , g(jv)|{0}), where |{0}) is the phonon vacuum. We consider two hnearly-polarized 
laser beams with wavevectors arranged such that ki — k2 is parallel to the molecular array. 
If the laser beams are far-detuned from any vibronic resonance the effective light-matter 
interaction operator can be written as V{t) — — 5'Ar[cJe~*'*'* -|- CqB^^], where qn is a size- 
dependent coupling energy proportional to the amplitudes of both laser beams, q — |ki — 
k2| and co = coi — lj2 are the net momentum and energy transferred from the fields to 
the molecules. For short interaction times (linear response), the system is excited from 
\g) into the one-particle sector with a probability proportional to the spectral function 
A{q,uj) — —lm.[G{q,uj)]/7r where G{q,uj) — {g\cq{huj — 1-1 + i'r))~^c^g\g) is the one-particle 
Green's function. As a result, for any q the polaron energy Eg equals the energy huu of the 
lowest-energy peak in A{q,uj), in analogy with ARPES measurements [12]. 

To measure A{q, u), the stimulated Raman excitation rate can be determined using state 
selective resonance enhanced multi-photon ionization (REMPl), see Fig. 3(b). With some 
probability this converts the rotational excitation (the "particle" ) into a molecular ion which 
can be extracted from the chain and detected by a multi-channel plate ion detector. Ioniza- 
tion and subsequent detection efficiencies for a 2 step REMPI processes can easily exceed 
20% [23] , and a properly gated integrator can resolve the arrival of a single molecular ion. 
Using a 3D lattice, a set of uncoupled parallel ID arrays can be realized and excited simul- 
taneously, increasing the signal to noise ratio of the detection step. 

In summary, we presented the first (to our knowledge) phase diagram for a mixed type 
polaron Hamitonian [30], which showed that polaron physics is much richer than previously 
thought and that sharp transitions may occur even for dominantly type-(i) Hamiltonians 
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FIG. 3: (Color online) (a) A two-photon stimulated Raman transition creates a polaron state with 
a well defined momentum q and energy co = coi — uj2- (b) The presence of the quasiparticle is 
subsequently detected using resonantly-enhanced multi-photon ionization (REMPI). 

(i? ^ 1). We showed that polar molecules trapped in optical lattices can be used to study 
this physics, and proposed an ARPES-like detection scheme to directly measure the polaron 
dispersion and thus identify the transitions expected to occur in such systems. 

Polarons with type (i) coupling have been observed for a single atomic impurity immersed 
in a Fermi gas [24-26], and proposed for realization in lattice setups using atom-molecule 
systems [27], self-assembled crystals in strong DC fields [18], and recently trapped ions [28]. 
Realization of type (ii) coupling has been considered using Rydberg atoms [29] . Using these 
systems it would be possible to explore the region of the phase diagram with kgs = for 
t < (or equivalently fcgs = tt for t > 0). Here we showed that using trapped molecules in 
weak DC electric fields it is possible to explore transitions into the phase with < |A;gs| < n, 
where novel polaron physics is expected to occur [8, 10]. 

Many other aspects of single polaron physics can be investigated with trapped polar 
molecules. Examples include investigating the effects of dispersive phonons (most theoreti- 
cal work assumes Einstein bosons), or novel effects resulting from quadratic particle-boson 
coupling in the strong coupling regime. Studying polaron phase diagrams in higher dimen- 
sions is easily achieved with the same experimental scheme. Generalizations to studies of 
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bi-polarons are also of significant interest, to understand the pairing mechanism for domi- 
nantly type-(ii) models (most bi-polaron studies are for type-(i) Holstein and Prohlich mod- 
els). Finally, one may also be able to adapt the polar molecules systems to study finite 
polaron concentrations and search for quantum phase transitions [31]. 

Acknowledgements: Work supported by NSERC and CIFAR. FH would also like to thank 
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Supplementary Material 



Rotational state dressing schemes 

Static field dressing - This dressing scheme involves an homogeneous DC electric field acting 
on a one-dimensional array of Af polar molecules. The field can have an arbitrary orientation 
with respect to the array axis, described by the polar angle ©. We consider molecules 
each prepared in its lower vibrational and electronic state. The rotational state of the 
molecules is described by the pcndular state \N,M), which is an eigenstate of and 

— BeN^ — dE-Dc cos ^, where N is the rotational angular momentum operator, is the 
rotational constant, d is magnitude of the permanent dipole moment, and 9 is the polar angle 
of the internuclear axis [? ]. We use the two-level subspace \g) — |0,0) and |e) = |i,0) to 
define the particle creation operator cj" = \ei){gi\. If the state |1, 0) is energetically separated 
from the degenerate states |1,±1) by an amount larger than the dipole-dipole interaction 
energy between neighbouring molecules in an array, the two-level approximation is valid and 
the dipole-dipole interaction operator reduces to Vdd = (1 ~ 3cos^ 6) do{i)do{j), where 

Tij is the intermolecular separation vector, © its polar angle with respect to the dc field, 
and dp is the p-component of the space-fixed dipole operator in spherical coordinates. This 
simplification of the dipole operator is accurate for electric field strengths £^dc > 1 V/cm, 
regardless of its orientation with respect to the molecular array. The dipolar energies ti2 = 
{e-igj\Vdd\giej) and D12 = {eig2\Vdd\eig2) - (5'i5'2|Vdd|fi'ifi'2) arc evaluated by diagonalizing 
numerically and using the eigenvectors to compute the matrix elements of Vdd for a given 
value of the field strength -Bdc- 

Combined static and infrared dressing - This schme involves a DC electric field as before 
and an additional infrared laser fields that induce Raman couplings between rotational states 
of the same parity. The Raman coupling scheme is illustrated in Fig. 1. It consists of two 
uncoupled lambda systems Ai = = 0)|0, 0), f = 1) |i, 0), I^E, f = 0)|2,0)} and 

A2 = {\^^,v^ 0)|i,0), I^E,^ = 1)|2,0), pE,^ = 0)|3,0)}. Here \^E,v) denotes the electronic 
and vibrational quantum numbers. Four linearly polarized infrared fields are needed to 
implement the scheme. The corresponding transitions are far detuned by Su > 10^ MHz 
from other rotational fines in the presence of a DC field Ej^c ~ 1 kV/cm. 

Let = \fl\e^'^ be the Rabi frequency associated with a dipole-allowed transition. For 
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FIG. 4: (a) Infrared dressing of rotational states in the presence of a dc electric field. Two 
independent lambda systems are used to construct the dark state superpositions \g) = ^/a\0,0) + 
71^12, 0) and |e) = Vb\i,0) + ^1 - 6|3,0). (b) Dipole-dipole energy ratio R = Di2/ii2 as a 
function of the DC fields strength -Bdc for DC field rotational dressing (a = 6 = 1) and combined 
AC-DC dressing with a = b= 1/2. 

each lambda system the adiabatic state with zero eigenvalue (dark state) can be written as 
[2] 

|Ao) = cosa;(t)|ei) — sina(t) exp(2/3)|e3), (3) 

where the mixing angle a{t) is defined by tanQ;(t) = \Qp(t)/Qsit)\, and the relative phase 
of the fields is P = (pp — (ps- The RWA holds when max {A, Qp, Qs} ^ B^. The two lambda 
systems in Fig. 1 are uncoupled due to the anharmonicity of the rotational spectrum for 
small detunings Ai and A2. Selection rules restrict the coupling between rovibrational states 
with different angular momentum projection M. 

We want to use the dark states of the lambda systems Ai and A2 to define the two- 
level subspace \g) and |e), respectively. It is in this field-dressed space where excitation 
hopping processes take place in a molecular array. Let the molecules be initially prepared 
in the rovibrational ground state \g) = \^T,,v = 0)|0,0). We assume that only the Stokes 
fields ^Siit) and ^Is-ii^) present. In this case the mixing angles and a2{t) vanish 

and the dark states are the same as the DC dressing scheme \g) = = 0)|0,0) and 

|e) = \^T,,v = 0)|1,0). By slowly turning on the pump fields fipi(t) and ^P2{t) the dark 
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states become 



1^) = cosai(t)|^E,'y = 0)|0,0) +sinai(t)|iE,'t; = 0)|2,0), 
|e) = C0SQ;2(t)|^S,'t; = 0)|i,0) +sinQ;2(i)|^S,'t; = 0)|3,0). 

The Raman lasers have the relative phases (3i — (32 — tt, which can be achieved experimen- 
tally. The states \g) and |e) can be manipulated independently by adiabatically tuning the 
pump laser intensities relative to the Stokes beams. We consider the specific situation where 
CKi — a2 — 7r/2. Omitting the labels v and M, the excitation hopping amplitude tij is then 
given by 

1 _ 3cos^ /l\ 

tij ~ j {dlo + + 4^ + 2d2idw + 2(^32^10 + 2^21(^32} (4) 

where dNN' — {N\do\N') — d^'N- We have ignored couplings with A7V > 1, which are 
suppressed in the presence of weak dc electric fields E-qc < -Be/ d. Equation 4 becomes exact 
in the limit i?DC where the transition dipole moments can be evaluated analytically 
using angular momentum algebra [l].The transition dipole moments in Eq. (4) have the same 
sign and therefore contribute to the enhancement of the hopping amplitude = f ti2 with 
respect to the pure dc dressing scheme. The diagonal energy D12 has a similar behaviour as 
in the type-I coupling only for weak dc electric fields dE-Qc/Be <^ 1. For non-perturbative 
fields, contributions from higher rotational states N >2 modify significantly its static field 
dependence. In Fig. 1(b) we show the achievable values of the ratio R = Di2/t\2 using 
the combined AC-DC dressing scheme described above, as function of the DC field strength. 
The hopping amplitude is enhanced by a factor / ~ 2 with respect to the static scheme. The 
ratio R has the same sign as the static dressing case {R < 0) for weak fields £^dc < ^B^/d, 
but is up to an order of magnitude smaller. The ratio changes sign at E'dc ~ 2.SBe/d and 
has a local maximum -Rmax ~ 0.32 at Ej^c — 7Be/d. 

Phonon dispersion.- In this work the dispersion of the phonon spectrum is ignored due 
to the current limitations of the Momentum Average (MA) approximation we use in the 
main text to solve the polaron problem [3]. For a ID array the phonon frequency can be 
written as Uq — UJo^/l + 12p^{q), where p — /Vo is the ratio between the ground state 
dipolar energy — {gi92\Vdd\gig2) and the optical lattice depth Vq, and 7(5') = 2sm'^{q) 
[4]. The dispersion is determined by the parameter p, which can be kept small for shallow 
lattices fa 18 — 20Eji in the static dressing scheme only for weak dc fields E'dc < Bg/d. 
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For higher fields we have |T42p^ > |^i2|, thus the phonon dispersion becomes an important 
energy scale. For combined AC-DC dressing, the phonon dispersion can be safely ignored 
since |V^2^| <^ |ti2| over the entire range of electric fields in Fig. 1. 
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